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Abstract 

By starting from the stochastic Hamiltonian of the three correlated spins and mod- 
eUng their frequency fluctuations as caused by dephasing noisy environments described 
by Ornstein-Uhlenbeck processes, we study the dynamics of quantum correlations, in- 
cluding entanglement and quantum discord. We prepared initially our open system with 
Greenberger-Horne-Zeilinger or W state and present the exact solutions for evolution dy- 
namics of entanglement and quantum discord between three spins under both Markovian 
and non-Markovian regime of this classical noise. By comparison the dynamics of en- 
tanglement with that of quantum discord we find that entanglement can be more robust 
than quantum discord against this noise. It is shown that by considering non-Markovian 
extensions the survival time of correlations prolong. 
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1 Introduction 

The unavoidable interaction of any realistic quantum system with its environment can destroy 
coherence between the states of a quantum system. This is decoherence, which can be intro- 
duced by various noisy models. Since, this phenomenon is deemed as one of the main obstacles 
to the realization of quantum information processing, recently, much attention has been paid 
on influence of it in quantum process [7, 8, 9]. Decoherence dynamics of entanglement and 
quantum discord, as two very useful resources to perform various quantum tasks, under the 
influence of environmental noises has been extensively discussed [1, 2, 3, 4, 5, 6, 13, 14, 15, 16]. 
When the environment correlation time is much shorter than the relaxation time of the system, 
i.e., environment has no memory, the Markovian approximation can be used to quantified the 
decoherence dynamics [10, 11, 12]. This approximation simplifies greatly the mathematics to 
solve the dynamics of open quantum system, but, in truth every environment is non-Markovian, 
i.e., with memory. The non-Markovian effect is a kind of dynamic backaction effect on the 
system due to pronounced memory effect of environment which could compensate the lost 
coherence of quantum system. Preservation of entanglement and discord can be caused by this 
characteristic of non-Markovian decoherence. 

Many investigator have lately used the well-known classical Ornstein-Uhlenbeck (OU) process 
as a model of non-Markovian decoherence [17, 18, 19]. The OU process has a long history 
in physics and plays a central role in the mathematical descriptions of Brownian motion and 
Johnson noise [20, 21]. In this letter, we present the exact solutions for evolution dynam- 
ics of entanglement and quantum discord between three spins, initially described by pure 
Greenberger-Horne-Zeihnger (GHZ) or W state, under the classical OU noise. As expected, 
the life of entanglement and quantum discord prolong due to the feedback dynamical of non- 
Markovian decoherence. The aim of this paper is to compare the robustness of entanglement 
with that of quantum discord against to the OU decoherence. We find that entanglement can 
be more strong than quantum discord under this noise. 

The structure of this article is as follows: First we describe the model of OU decoherence and 
two correlation measures which we used for this study. In next section we present our results 
for entanglement and quantum discord between three spins with initial GHZ and W state. 
Summary of this letter are given in the last section. 
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2 The Model 

2.1 The Ornstein-Uhlenbeck decoherence 

The random fluctuations of a spin-frequency which are sufficiently small, can approximately 
be described as an OU process [17, 22]. Based on this, the stochastic Hamiltonian of the three 
correlated spins with random fluctuations frequency 

Htot{t)^j:lni{t)a^;\ (2-1) 

i=l ^ 

can be assumed as the Hamiltonian of interaction between system and the local OU noise. 
Under this condition, which is the independent transition frequency of the i-th spin, 

describe the local OU noise with the statistical mean value properties as 

M{^,{t)\ = 0, 

M[Q,(s)Q,(t)] = ^(s - = ^r,7e-^l^-*l . (2-2) 

In this classical noise, 7 is the bandwidth of the noise, which is related to the correlation time 
as Tc = ^. Figuratively, the correlation time can be deemed as the time scale over which the 
environment has memory and after it the environment is back to equilibrium and effectively 
has no memory. The coupling strength between the i-th spin and the environment is given by 
Fj, (i = 1, 2, 3), which for simplicity we assumed Fj = F for all spins. 
The time evolution of the total system can be calculated explicitly from 

Ptotit) = Uit)ptotiO)U\t), (2-3) 

where the stochastic unitary operator U (t) is the explicit solution for the stochastic Schrodinger 
equation as 

U{t) ^e-'Io"*"*^'^'^'. (2-4) 

Clearly, the unitary operator U (t) is dependent on the noise. By taking the ensemble average 
over the noise field, ( i.e., from the statistical mean ) the reduced density matrix of open spins 
system ps{t) is obtained 

Ps{t) ^ M[ptot{t)]. (2-5) 

So, the dynamics of the system density matrix can be described in terms of Kraus representa- 
tion as 

p,(t) = ^X,(t)p,(0)X/(t). (2-6) 



4 



By using the fact that Tr[ps{t)] = 1, a condition on the Kraus operators can be obtained as 
J2iKj(t)Ki(t) = I. In our case, the Kraus operators describing the interaction with the local 
OU noise are given by 



Ki = Fi (g) Fi, K2 = F^0Fi®F2, i^g = i^i ® ® i^i , = F2 ® Fi (g) Fi, 

Kr, = Fi® F2® F2, Kq = F2®Fi®F2, Kt = F2®F2®Fi, Kf^ = F2 ® F2 ® F2, 



where 



(2-7) 



-(T :)• -(f 

The parameters appearing in Fi and F2 operators are 

liit) = exp[f P{s-t)d's]=exp[^{t + -{e-^'-l)}], 
Jo z 7 

u{t) = ^l-fi'it)- (2-9) 

For sufficiently large values of bandwidth, 7 — )■ 00, i.e., the correlation time Tc — t- 0, we 
get j3{s — t) = r6{s — t), and hence the Markovian dynamic of OU noise is recovered. For 
this limit, the factor of decoherence is reduced to fi{t) exp[^^], hence, the coherence 
decay rate determine by the coupling strength between the spins and environment F. If the 
dynamics of system and environment is such that the correlation function in Eq. (2) can not 
be replaced by a delta function, (i.e., 7 — >■ ), the dynamics is non-Markovian and memory 
effects of the environment have important roles. In this case, by using the approximation 
Q-jt I _ ^-f-j^ 1(7^)^, the factor of decoherence, can be expanded as /x(i) — >■ exp[^7rt^]. 



2.2 Measuring entanglement 

Since entanglement is conceived as a resource to perform various tasks of quantum information 
processing [23, 24, 25, 26], knowledge about the amount of entanglement in a quantum state 
is so important. Indeed, awareness from the value of entanglement, means knowing how well 
a certain task can be accomplished. The quantification problem of entanglement only for 
bipartite systems in pure states [27] and two-qubit system in mixed state [28] is essentially 
solved. In multi-partite systems, even the pure state case, this problem is not exactly solved 
and just lower bounds for the entanglement have been proposed [29, 30, 31, 32]. Here, in order 
to determination the exact minimum of entanglement between three spins, we use the lower 
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bound of concurrence for three-qubit state which is recently suggested by Li et al. [33] 

rsip) = ^AUcrr + {Cfr + iCfy)l, (2-10) 
v3 j=i 

12|3 

where Cj is terms of the bipartite concurrences for qubits 12 and 3 which is given by 

Cf^ = max{0, Af l'(l) - Af ''(2) - Af ''(3) - Af ''(4)}. (2-11) 

12 13 

In this notation, A^ (k), {k = I. A), are the square nonzero roots, in decreasing order, of 
the non-Hermitian matrix ppj^'^. The matrix pj^'^ are obtained from rotated the complex 
conjugate of density operator, p*, by the operator S'j^'^ as p]^'^ = S'J^'^ p* S'J^'^. The rotation 
operators S'j^'^ are given by tensor product of the six generators of the group S0(4), {Lj'^), and 
the single generator of the group S0(2), (Lq) that is S'J^'^ = Lj^ L^. Since the matrix S'j^'^ 
has four rows and columns which are identically zero, so the rank of non-Hermitian matrix 
pp]^'^ can not be larger than 4, i.e., A]^'^(k) = for k > 5. The bipartite concurrences C^^'^ 
and C^^l-*^ are defined in a similar way to C^'^^^. 



2.3 Measuring quantum discord 

Another kind of quantum correlation which is in general different from entanglement has been 

designated as quantum discord [34]. This measure of quantum correlation, which is arising from 

the difference between two quantum extensions of the classical mutual information, proved its 

abilities as a fundamental resource for quantum information tasks [35, 36, 37, 38]. Recently, 

many efforts to generalization of quantum discord to multi-partite systems have been made 

by different authors [39, 40]. The global measure of quantum discord which is obtained by 

a systematic extension of the bipartite quantum discord is the result of such efforts [41]. By 

using this measure one can quantifies the quantum discord of an arbitrary multi-partite state. 

Under a set of von-Neumann measurements as 

cos^d) e^^'cos(|)sin(|)\ 

^ Ve~'^'cos(|)sin(|) sin^d) )' 



n«^( . , -e-cos(|)sin(|)- 

2 I _emcos(|)sin(|) cos^d) ^ ^ 



which rotated the direction of the basis vector of Z-th spin with 9i e [0, tt) and (pi e [0, 27r), the 
global quantum discord has the form 

Dip) = min[5(p||$(p))-i:5(p(ol|$«(pO)] 
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= min [5($(p)) - S{p) - - S{pi))]. (2-13) 

In this notation, S{p) = —Tr[p\og2 p] is the von-Neumann entropy of the density matrix and 
$(')(p;) = nf-'-l-n^'^pi n^'-* is the reduce density matrix after performing the measurement 

on the Z-th spin. The matrix obtain by carrying out the 8 projective measurements 

n4 = n^^^ ® nf ^ ® nS^^ , = u'^'> ® n^^^ ® , ng = n^^^ ® nf ^ o ^ , 
Uj = n^^^ ® n^^) nf \ ng^n^^^^n^^^^nf , (2-i4) 

on the three correlated spins as: $(p) = Sm=i P ^m- In order to ehminate the dependence 
of quantum discord on the measurement operators, we must find the measurement basis that 
minimizes D{p). 

In the following wc will investigated the evolution of entanglement and quantum discord, as 
two different kinds of the quantum correlation, between three correlated spin which initially 
described by pure GHZ or W state under the OU Markovian and non-Markovian decoherences. 

3 The dynamics of three correlated spins under the OU 
noise 

In this section we present time evolution of entanglement and discord between three spins 
which is coupled with OU classical noise. We assume the open system is initially prepared 
in inequivalent class of pure three qubit state with maximally quantum correlation which is 
known as GHZ and W states. 

3.1 Initial GHZ state 

First, we suppose that the initial state of the system is pure GHZ state 

\GHZ) = ^(|000) + |111)). (3-15) 
v2 

According to Eq. (6), the reduced density matrix of the system under the OU decoherence 
can be expressed as 

PGHzit) = ^{|000)(000| + |111)(111|} + ^{|000)(111| + |111)(000|}. 
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(3-16) 

The concurrence of the above density matrix can be easily computed as 

n(pGHz(t)) = l^'itMpGHziO))- (3-17) 

Now, we turn to calculating the quantum discord from Eq. (13). By tracing out any two spins 
of the density matrix (16), one can obtain the reduced density matrices of the subsystems 
as Pi = P2 = P3 = I- So, Y4=i{S{^^''\pi)) - S{pi)) = 0. The von-Neumann entropy of the 
PGHzit) can be computed as 

S{pGHz{t)) = 1 - ^i|^log2(l + p^(i)) - ^^|^log2(l - (3-18) 

After some calculation to find the desired measurement basis which minimizes the quantum 
discord, we deduced that the eigenbasis of the Pauli matrix are the best measurement basis. 
Under this local measurement we obtain S{^{pGHz{t))) = 1 and 

D{pGHz{t)) = ^^^|^log2(l + p-^(0) + log,(l-/x-^(t)). (3-19) 

In order to comparison the robustness of the entanglement (17) and the quantum discord 
(19) against to OU decoherence, we plot these quantities as a function of the dimensionless 
scale rt for both Markovian (a) and non-Markovian (b) regime in Fig. 1. It is explicit 
that, entanglement can be survived under the Markovian and non-Markovian conditions more 
than the quantum discord. Moreover, observe that due to the influence of memory effect of 
environment the entanglement and the quantum discord have longer life. In the other words, 
one could imagine that non-Markovian effect is a key factor to preserve quantum correlation 
since non-Markovian effect is a kind of backaction effect which could compensate the lost 
coherence of quantum system. 

3.2 Initial W state 



When the three correlated spins is initially prepared in the W state 

\W) = ^(|100) + |010) + V2|001)), 



(3-20) 
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the density matrix dynamics, according Eq. (6), can be express as 



/O 




Pw{t) = 




2 



0\ 
yriy?{t) y/2[i\t) 



V2//2(i) 


^^2il\t) 


VO 



1 



















1 

















0/ 



(3-21) 



From the lower bound of concurrence for three-qubit state Eq. (10) for the above matrix we 
obtain 



(3-22) 



In order to determine quantum discord according to Eq. (13), we first evaluate the one-spin 
density matrices representing the individual subsystems by tracing out two spins. These reduce 
density matrices have the flowing form p-^ — p2 — p-i — |{3|0)(0| -|- |1)(1|}. Next, by using 
the same procedure done in previous section to find the desired measurement basis which 
minimizes the quantum discord, we obtain the minimum of quantum discord with adopt local 
measurements in the eigenbasis for any spin. Under such measurements the state of the 
single spin is not induced therefor we have Ya=\{^{^^^^{Pi))~ ^{Pi)) = S{^{pw{t))) = §■ 
Finally, by inserting the von-Neumann entropy of the time-dependent reduced density matrix 
for three correlated spins in Eq. (13), we obtain the flowing form for quantum discord 

^{Pw{t)) = ^{h -r ll\t)) -r\{\ - p\t))\0^^{\ - p\t)) 



log2(3 + p'it) - ^1 - 2p\t) + np \t)) 
+ (3 + /i2(t) + v/l - 2/i2(t) + \lii\t)) 



log2(3 + p'{t) + Jl- 2p?{t) + nii\t))}. 



(3-23) 



Notice that we shall re-normalize the Eq. (23) in such a way that we have D{pw{^)) — 1- 
From the expressions of entanglement (22) and quantum discord (23), it can be found that 
for pit) quantum correlation between three spins disappear. In Fig. 2, the dynamics 
of these quantities for both Markovian (a) and non-Markovian (b) dynamics are illustrated. 
Clearly, robustness of entanglement versus OU Markovian and non-Markovian decoherence is 
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more than quantum discord. 

Notice that, since the factor of decoherence fi{t) appears as quadratic in quantum correlation 
of W state, this state has hfc longer than GHZ state with cubic decoherence. Comparison of 
plots in Fig. 1 and Fig. 2 justifies this fact. 

4 Summary 

In summary, we have investigated the exact decoherence dynamics of quantum correlations, in- 
cluding entanglement and quantum discord, between three correlated spin with initial GHZ and 
W state under a local dissipative OU process. By studying both Markovian and non-Markovian 
regime of this classical noise, we have found that under Markovian dynamics quantum cor- 
relation suffers sudden death, while due to the influence of memory effect of environment in 
non-Markovian dynamics, it has longer life. Moreover, comparison of the survival time of 
entanglement and quantum discord make it clear that against dephasing OU decoherence, en- 
tanglement is more robust than quantum discord, because it decays exponentially while the 
discord of the same state decays logarithmic. Also, our results have shown that W sate has 
strong dynamics under decoherence than GHZ state. 
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Figure 1: (Color online) Entanglement (solid blue) and quantum discord (dashed red) dynamics 
of three correlated spin with initial GHZ stat under OU decoherence. Here, we have chosen 
the parameters F = 1. Fig. (a) corresponds to Markovian regime with fi{t) exp[^^] and 
Fig. (b) to non-Markovian regime with ii{t) exp[^7rt2] and 7 = 0.01. 
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Figure 2: (Color online) Entanglement (solid blue) and quantum discord (dashed red) dynamics 
of three correlated spin with initial W stat under OU decoherence. Here, we have chosen the 
parameters F = 1. Fig. (a) corresponds to Markovian regime with /x(t) exp[^|^] and Fig. 
(b) to non-Markovian regime with ii{t) exp[^7Ft2] and 7 = 0.01. 



